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Abstract 

We discuss questions pertaining to the definition of 'momentum', 'momentum space', 'phase 
space', and 'Wigner distributions'; for finite dimensional quantum systems. For such systems, 
where traditional concepts of 'momenta' established for continuum situations offer little help, we 
propose a physically reasonable and mathematically tangible definition and use it for the purpose 
of setting up Wigner distributions in a purely algebraic manner. It is found that the point of view 
adopted here is limited to odd dimensional systems only. The mathematical reasons which force 
this situation are examined in detail. 
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1. INTRODUCTION 



The description of states of quantum systems with the help of Wigner distributions, 
in the case of continuous Cartesian position and momentum variables , has found 

applications in a wide variety of contexts. These include the original aim of calculating 
quantum corrections to classical statistical mechanics, various aspects of quantum optics, 
and quantum chemistry Q]. 

There has been for some time considerable interest in developing a similar method 



12 1 . In the works of 



for quantum systems possessing a finite dimensional state space 
Feynman|6| and of Woottersjsl, which constitute some of the early efforts in this direction, 
it was found necessary to treat independently the cases of two-dimentional state space, and 
of an odd prime dimensional state space. In the former, corresponding to the case of a 
qubit, the insights from the treatment of a spin 1/2 angular momentum in quantum me- 
chanics were used effectively. The treatment of the odd prime dimension case has been 
geometric in spirit. From the start it is assumed that the phase space is a square array of 
points, namely that 'momentum space' is of the same 'size' as 'position' or configuration 
space. In this phase space, the notions of straight lines, parallel lines, and foliations or 
striations are set up and used in stating the properties desired of a Wigner distribution and 
then a solution is proposed The arithmetic of a finite field of odd prime order is 

used in the constructions. For general odd dimensions or a power of two, one has to take 
the Cartesian product of the elementary solutions. 

The purpose of the present work is to provide an alternative treatment of this problem 
based on an algebraic approach. As we will find, apart from reproducing known results (in 
the odd dimension cases), essentially new possibilities will also come to light. Our work is 
motivated by the following points of view. In elementary mechanics, the primitive notion 
of momentum is as a numerical measure of rate of change of position in space, the factor 
of proportionality being the mass. In canonical mechanics we move to an algebraic role 
for momentum, as generator of translations in configuration space, with mass not playing 
an explicit role. In Cartesian quantum mechanics this algebraic role is expressed by the 
canonical Heisenberg commutation relations and is fundamental. 

Returning to classical canonical mechanics usually the configuration space Q is taken to 
be a differentiable manifold of some dimension. Then the concepts of momenta and phase 
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space are given immediately by the automatic construction of the cotangent bundle T*Q. 
There are of course special features in the Cartesian case Q = M", not shared by the general 
case. Then T*Q = R^", and both positions and conjugate momenta are elements of real 
linear vector spaces of the common dimension, n. Thus they are of the same algebraic nature, 
allowing us to form real linear combinations of them in a physically significant manner. This 
leads to an important role for the groups Sp{2n, R) in the classical case, and to the double 
covers Mp{2n) after quantisation. For a general configuration space Q different from M", 
for instance or the manifold of a Lie group one sees that coordinates and momenta 
have intrinsically different algebraic natures. Locally momenta always belong to a linear 
vector space, while positions may have no such characterisation. In these cases, of course, 
quantum mechanics uses suitably square integrable functions on Q as wave functions, and the 
basic operator relations are inspired by what one knows about the classical T*Q. However, 
already in the simple case of Q = S^, a particle confined to the surface of the unit sphere, one 
sees an interesting new possibility emerging. Since T*§^ is not parallelisable, the canonical 
momenta cannot be globally defined as variables. However one can view §^ as the coset space 
S0{3)/ S0{2), and so think of canonical momenta as the constrained components of a three 
dimensional angular momentum, obeying the algebraic requirement of being orthogonal to 
the (unit) position vector. This makes the momenta globally well defined, besides bringing 
into the picture a role for the group 5*0(3) acting transitively on 

We now turn to an A^-level quantum system. The state space is a complex A^-dimensional 
Hilbert space, wave functions are complex A^-component column vectors, but properly speak- 
ing no notion of canonical momenta is given in advance in any intrinsic sense. There are 
no fundamental operator relations analogous to the canonical Heisenberg commutation re- 
lations. In this case, classically we may view the configuration space Q as a finite set of 
discrete points, without any concept of canonical momenta. This leads to the viewpoint: for 
finite level systems, the problems of defining momenta, momentum space and phase space 
should be regarded as part of the problem of setting up Wigner distributions in a physically 
reasonable manner. 

From this discussion we extract and explore the idea that for any choice of configuration 
space Q, whether a differentiable manifold or a finite discrete set, the concept of momentum 
is any group G that acts transitively on Q. In case Q is a differentiable manifold and G acts 
on Q via diffeomorphisms, we move a step closer to the usual meaning of momentum by 
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consideration of the generators of these diffeomorphisms subject to as many hnear algebraic 
constraints as the difference in the dimensions of G and Q. 

This general viewpoint means that Q is a coset space G/ H with respect to some subgroup 
H G G. We will however assume for definiteness that H is trivial, i.e., Q = G. Thus the 
configuration space is itself a group. We will find that using the representation theory of G, 
in the quantum case this leads to an algebraic way of defining momenta and phase space. 
However when we turn to the problem of defining Wigner distributions the method is limited 
to odd dimensions only for a reason that will emerge later. 

This work is organised as follows. In Section 2, following the line of thought outlined 
above, we take a group G of order N as the 'classical' configuration space underlying an 
A^-level quantum system; and proceed to set up the necessary machinery for introducing the 
notions of 'positon space', 'momentum space' and 'a point in phase space', recapitulating, 
in the process, relevant aspects of the representation theory of finite groups. A noteworthy 
feature of the phase space thus obtained is that it is always an x array for any A^ 
regardless of the nature of G. Having developed the notion of 'phase space' for an A^-level 
quantum system, for any A^, in Section 3, we take up the question of constructing a Wigner 
distribution on such a phase space. This requires introducing the concept of a 'square root' 
of a group element, which, as it turns out, is possible only when A^ is odd. Restricting 
ourselves to such dimensions, we show that, in general, there are two equally good ways of 
defining Wigner distributions which coincide only when G is an abelian group. Choosing 
one of them, and keeping G general, we construct the corresponding phase point operators 
j^,^^ and find, in particular, that the phase point operator at the origin, like its continuum 
counterpart, has the satisfying feature of being an inversion operator. We then discuss in 
some detail the special case when G is abelian and analyse the two possibilities (a) A^ an 
odd prime (b) A^ an odd composite. While in the former case, G gets uniquely fixed to 
be a cyclic group of order A^, in the latter case, in general, one has more than one choice 
for G leading to distinct definitions of the Wigner distribution. (The material presented in 
this section makes use of the results in though with some changes of notation). In 

Section 4 we turn to the case when G is nonabelian and analyse the relation between the 
two choices for the Wigner distribution. We further show that it is possible to define an 
'extended' Wigner distribution which contains the two Wigner distributions as its marginals. 
However, as is explicitly demonstrated, this can only be achieved at the cost of introducing 
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overcompleteness. Section 5 contains our concluding remarks and further outlook. 



2. FINITE ORDER CONFIGURATION AND PHASE SPACES 

Our aim is to describe general states of an A'"-level quantum systems (for odd A^) via 
Wigner distributions. The space of (pure) quantum states is an A^-dimensional complex 
Hilbert space T-C^^\ each vector \ip >G Ti.^^^ being an iV-componcnt complex column vector. 
The corresponding or comparison classical configuration space Q is a finite discrete set of N 
points. We wish to regard Q as a finite group G of order N. This makes it natural to label 
the points of Q as gi, g', - ■ ■ , a, b, ■ ■ ■ & G, with a distinguished 'origin' corresponding to 
the identity element e e G. In turn in the quantum case we introduce an orthonormal basis 
of 'position eigenstates' \g > for Ti^-^^: 

< g'lg >= Sg'^g, g, g' e G; 

Y^\9><9\=\. (2.1) 

A general l^/' >G HS^^ has a position space wave function 

i){g) =< g\i) >, 

<i,\i,>=\\i,\f = Y,\i,{g)\\ (2.2) 

9 

To move onto the concepts of momentum and phase space, we consider the transitive action 
of G on Q = G, namely of G on itself. Keeping in mind the general non-abelian situation, 
we have two such actions, namely left translations Lg and right translations Rg. 

{Lg,i,){g) = i^ig'-'g), 

Lg'\g > = \g' g >; 
{Rg'^)(g) = ^l^igg'), 
Rg'\g> = \gg'^^>- (2.3) 

(In the abelian case, Lg — Rg). Both Lg and Rg are unitary, and they obey 

Lgl Lg = Lgig, 

L^g' L^9 ~ Ltgig, 

Lgl Rg = Rg Lgl. (2.4) 
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Thus we have two commuting representations-the regular representations of G on Ti^^). We 
arrive at the notion of 'momentum eigen states' by performing their simultaneous complete 
reduction into irreducibles. For this we recall well known facts regarding the inequivalent 
unitary irreducible representations (UIR's) of G. We label them with a symbol j, the 
corresponding dimension being Nj. With respect to some chosen orthonormal basis within 
each UIR, we have corresponding unitary representation matrices 

D^{g) = {DlM), m,n=l,2,.--,iV,. (2.5) 

The following well known properties and facts will be needed: 

Composition Y^^'mwU) = DUg'g); (2.6a) 

n' 

Unitarity D^M* = (2.6b) 

n 

Orthogonality ^Dj^^{g) D^^,^,{g)* ^ ^ 5jj, Smm' Snn'', (2.6c) 
Completeness ^ TV . D^^ig) D^gT = ^1^'') 

jmn jm 

= Nd,,,r, (2.6d) 

J2^j = ^- (2-6e) 

j 

It is now natural to define an orthonormal basis of 'momentum eigenstates' for 7i^^\ com- 
plementary to the basis {\g >}, as vectors labelled \jmn >: 

\jmn> = )J^J2^^rnni9)\9>, 

9 

<g\jmn> = ^ ^ Dl^(g). (2.7) 



They obey 



< j'm'n'\jmn > = 5m'm ^n'n, 

\jmn < jmn\ = I, (2.^ 



jmn 

while the completely reduced actions of Lg and Rg are : 



Lg \jmn > ^mm'{9 ^) \jm'n >, 

m' 

Rg \jmn > = Y1 K'n{9) \jmn' > . (2.9) 



Therefore 'momentum' is 'conjugate' to 'position' in the sense that when jmn and g are 
given, a complex number ^jNj/N Dl^^{g) is determined. 

A vector \ip >E Ti^^^ thus has two complementary wave functions: ip{g) of eq. (j2.2p . and 
a momentum space wave function 

iJjmn =< jmn\ijj > . (2.10) 
If I?/' > is normalised, then it determines two complementary probability distributions, as 

m\' = Yl \^(9)\' = J2 \^^n^n\' = l. (2.11) 
g jmn 

For purposes of setting up Wigner distributions, we thus regard a 'point in phase space' as 
a pair {g;jmn). By eq. ()2.6e|l we then find that the 'size' of momentum space is precisely N, 
the same as the 'size' of Q: so the points {g; jmn) of phase space form a square array. We 
obtain this result here not as a matter of definition but as a consequence of an important 
result of group representation theory. In the phase space there is a distinguished 'origin' 
corresponding to (e;000): here g = e, and j = m = n = denote the trivial or identity 
representation of G. 

The developments of this section so far are valid for all A^. Further, each way of regarding 
Q as a group of order leads to a corresponding way of building up momenta and phase 
space. Thus if G and G" are two non isomorphic groups of order A^, they lead to two distinct 
notions of momenta and phase space. 

3. WIGNER DISTRIBUTIONS 

We begin by recalling an important property available only for odd A^, i.e., a property 
of finite groups of odd order. It is that in such a group each element g E G has a unique 
'square root' element written as or g^f"^: 

geG^ ^eG : {^f = g. (3.1) 
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The uniqueness and other useful properties of this notion are as follows: 



91 = 92 = ^/92, 9i = 92^ (3.2a) 



(v^)-' = ^/F'\ (3.2b) 
9^ = V99; (3.2c) 



9-'V9 = vV^; (3.2d) 
V99'9-' = 9\^9-'- (3.2e) 



Therefore in any sum over G of a function f{g) belonging to a linear space we have the 
freedom to write the sum in different ways: 

E/(^) = E/(v^) = E/(^')- (3-3) 

9 9 

A further consequence is : 

a = 9'~^9, b = g'g 
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g' = Vba ^, g = V ab = Vba ^ a; (3.4a) 
a = 99'^, b = gg <^ 

g' = Va-^b, g = bVb-^a = aVa'^b. (3.4b) 

In both cases, for given a and b, g and g' are unique. In the nonabelian case, these two 
equations differ, while in the abelian case they coincide. We will find that these two prop- 
erties of any group of odd order are essential to be able to define Wigner distributions with 
desirable properties using the present method. 

We have seen in Section 2 that a point in phase space is a pair {g; jmn). It is then natural 
to look for a Wigner distribution in the form W^{g] jmn) for given normalised >G Ti^^'^ 
and require that it have desirable properties with regard to reality, transformations under 
left and right group actions, reproduction of the two marginal probability distributions in 
eq. ()2.11|) . and traciality. Keeping in mind the general case of nonabelian G, it turns out 
that there are two distinct but equally good ways of setting up Wigner distributions, which 
reproduce complementary features of the momentum space probability distribution iV'jmnp. 
We define them in sequence. 

Wigner distribution I 
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For any given normalised l^/" >G Ti^^'^ we define a corresponding Wigner distribution by 
W^{g-3mm') = J2H9'''g)Dl„,ig' 'Mg'g)*. (3.5) 

9' 

This definition extends immediately to a general state p as 

W^ig;jmm') = J2< g'-'gWg > D^^'ig' (3.6) 

9' 

however for simplicity we shall work mainly with the pure state case. Using the properties 
()2.6I3.3|3.4|) as necessary, we find : 

W^{g;jmm'y = W^{g]jm'm); (3.7a) 
^' = Lg,ilj : W^-{g;jmm') = ^ Dl^„^^{g')W^{g'-^g;jmim[) x 

mi ,m'j^ 

D' , ,(g'-^), 

< = Rg>^ : W^.{g;jmm') = W^gg'; jmm'); (3.7b) 

9 " 

^^iV, l^4^?;jmm) = \^Pig)\'; (3.7c) 

jr'm 

^ ^ ^ iV, iy^((7; jm'm) W^g^jmm') = \ < (f)\^ > (3.7d) 

g jmm' 

From the last traciality property it is clear that {ip > is completely described by W^{g] jmm') 
as far as all quantum mechanical probabilities are concerned. For the marginals, we see that 
the 'position space' probability distribution |^((7)P is faithfully reproduced. However, in 
the nonabelian case, the 'momentum space' probability distribution (in a polarised form) 
is reproduced in what we may call a 'coarse grained' version. This is further discussed in 
Section 4. 

Properties ()3.7al3.7bp and the second of ()3.7cp are evident essentially upon inspection. 
The first of eq. ()3.7c|) follows upon use of eq. ()3.4ap and changing summation variables. For the 
traciality property ()3.7d)l slightly more work is needed. Using eq. ()2.6d)l at the appropriate 
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stage, and then eq. ()3.4ap . we get: 

^ NjW^{g; jm'm)W^{g; jmm') = 

g jmm' 
g,g' ,g" jmm' 

= J2 <Pi9"'9m9"9m9"''9)<l>i9'9r^",g'-^ (3.8) 



9,9 ,9 



which immediately gives ()3.7d|) . 

Wigner distribution II 

Now in place of eg. ()3.5p we define a different (in the nonabelian case !) Wigner distribu- 
tion by 

W^{g-,3nn') = J2^igg'-')DUg' 'mgg'T- (3-9) 

9' 

Then in place of eqs. fl3.7p we find : 

W^{g-jnn'y = W^'{g; jn'n); (3.10a) 
= Lg,i) : W^'{g;jnn') = W^'{g'~^g; jnn), 
r = Rg'^:W^,.{g-jnn') = DL^(<7')H^/(^?^?'; jnin;)D^,„,(<7'-i); 

(3.10b) 

9 rn 

lj2N,W^'i9;jnn) = mg)f; (3.10c) 



N 

jn 

J^J2J2^^K(9;jnn) W^'{g;jnn') = | < > \'. (3.10d) 

g jnn' 

Having exhibited the existence of these two ways of defining the Wigner distribution 
distinct for nonabelian G, we hereafter mainly work with definition I, except in Section 4 
where we explore the relationships between the two. 

Phase Point Operators 

In analogy with earlier treatments, it is possible to express the Wigner distribution (j3.5p 
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as the expectation value of a 'phase point operator' in the state IV" >: 

W^{g; jmm') = < ilj\W{g; jmm')\ilj >, 

W{g-3mm') = \g'g > Dl^,{g' < g'-'g\. (3.11) 



Since with the help of eq. ()2.6d|) and ()3.2a|) we are able to 'invert' this to read 

N 



g'g X g'-'gl = E § ^'mm'ig' 'rW{g;jmm'), (3.12) 



jmm' 

we see that they provide a basis for the space of all operators on T-C^^\ They have the 
specific properties: 



W{g; jmm'Y = W{g] jm'm); (3.13a) 

N ^ 



-j^ W{g; jmm') = '^^\jmn >< jm'n\; (3.13b) 



g n 

J^J^'^^S-Jmm) = I^X^I; (3.13c) 

jm 

Ti[Wig;jmm')] = 6^^r, (3.13d) 
TT[W{g;jmn) W{g'; j'm'n')^ = —6g^g'6jf6mm'Snn'; (3.13e) 

which are easily obtained by eqs. ()3.4al 2.71 2.6(^1 as appropriate. 

We noted earlier that the point (e; 000) in phase space may be regarded as a distinguished 
'origin'. The corresponding phase point operator W{e; 000), from eq. ()3.1H) . has the simple 
form and action 

l?(e;000) = J2\9'><g"\ 

g' 

l?(e; 000)1(7 > = 1^?"' > • (3.14) 

We may call this as the inversion operation. This is the analogue, in the present finite 
dimensional case treated algebraically, of the known result in the case of continuous variables 
that the phase point operator at the origin of phase space is a multiple of the parity operator. 

The Odd Prime Case 

In case is an odd prime, the only choice for G is to be an abelian cyclic group of order 
iV, so we exhibit it as 

G = {e,g,g',---g''-V = e} 

= {<7ln = 0,l,--- ,iV-l}. (3.15) 
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The UIR's are all one dimensional labelled by an integer j = 0, 1, 2, ■ ■ ■ , — 1: 

The unique square root of g"' is easily computed: 

^ = gin+N)/2^ n = l,3,5,--- ,iV-2; 

= n = 2,4,6,---,iV-l. (3.17) 

The two Wigner distribution definitions ()3.5I3.9|) coincide: 

Af-l 

W4g^;j) = J2 ^(^/'^-">'^^^-'"'/^^(^"+"')*, (3.18) 

n'=0 

reproducing a known result P]. There is no 'coarse graining' involved in the recovery ()3.7c|) 
of the marginal momentum space probability distribution. 

For general odd A^, we may have again a unique G in the abelian case, or more than one 
distinct abelian choice. For example, for = 9, the two choices G = Cg and G" = C3 x C3, 
where C„ is the cyclic group of order n, are nonisomorphic, so they lead to distinct definitions 
of the Wigner distribution. On the other hand for = 15, C15 and C3 x C5 are isomorphic. 
The general pattern can be seen along these lines. 

The first case of a nonabelian G occurs when A^ = 21, some others occur for instance 
when A^ = 27, 39, 45, 55, 57, 63, 75, ■■■ . It is in these cases that the 'coarse graining' of the 
momentum space probability distribution seen in eqs. ()3.7cf 3.1()c|l in two different ways, is 
involved. We examine this in some detail in the next section. 



4. ANALYSIS OF THE NONABELIAN CASE 

Let A^ be an odd nonprime such that a nonabelian group G of order A^ exists. In that 
case some further interesting properties of the two definitions of Wigner distributions given 
in Section 3 can be exhibited. To begin with, by a simple change of variable g' g^^g'g in 
equation eq. ()3.9|) we can put the second definition into the form 

W^\g;jnn') = Y^i^ig'^' 9)Kn'{9~' 9' '9)^9 9)*. (4.1) 
3' 
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Comparing with W^{g\ jmm') in eq. ()3.5p . and for each g and j treating the two Wigner 
distributions as Nj x Nj matrices, we see that 

W^{g-3) = D^ig-')W4g;j)D^ig), (4.2) 

making clear the relation between them in detail. 

Next we turn to the coarse graining of the momentum space probability distribution 
involved in eqs. ()3.7cl3.10cp . It is in fact possible to define an overcomplete or redundant 
Wigner distribution in the following manner: 

w^{g;jmn'm'n) = '^'^{g'~^g)Dl^^,{g'g)Dl^^,{g-^g')'^{g'g)* 

9' 

g' 

From here the two actual Wigner distributions emerge by partial tracing: 

W^{g; jmm') = (g; jmn m'n) , 

n 

W^'{g;jnn') = ^w^{g;jmn mn). (4.4) 

m 

The quantities w^{g]imn' m'n) have the properties 

w^{g; jmn' m'n)* = w^{g; jm'n mn'); 
-j^2^w^{g]jmn' m'n) = ^pjm'ni^jmn''^ 

9 

j^^NjW^{g;jmnmn) = \ip{g)\'^. (4.5) 

jmn 

These results are consistent with eqs. fl4.4l 3.71 3.10|) . and the behaviours under left and 
right group actions can be easily worked out. We now see that all the momentum space 
probabilities (in a polarised form), without any coarse graining, can be recovered from the 
distribution w^{g] jmn' m'n). However this distribution is overcomplete, as is demonstrated 
by showing explicitly that it can be reconstructed from, say, W^{g; jmm'). Combining the 
definition ()4.3p with eqs. ()3.11l3.12|) we have: 

w^{g;jmn' m'n) = ^ < -^{g'g >< g'-^glip > D^^^,{g'g)Dl^,{g-^g') 

9' 

g' j"m"n" 

W4g-j"m"n"). (4.6) 
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Again it is easily checked that this is consistent with eqs. ()4.2I4.4|) . This clarifies the 
situation in the nonabelian case. The full momentum space probability distribution \ipjmn\'^, 
more generally the polarised expressions 4'jm'n4'jmn'y obtainable in principle from the 

Wigner distribution Wjp{g] jmm'), even though at first sight eq. (j3.7c|) seems to suggest that 
only the coarse grained expressions i'jm'n^'jmn '^^^ be recovered. It is indeed the case 
that W^{g; jmm') describes \ip >< i'l completely. 

5. CONCLUSIONS 

In this work we have proposed a way of associating a phase space with an A^-level quantum 
system by taking the corresponding 'classical' or comparison configuration space as consist- 
ing of the N elements of a group G of order A^. A gratifying feature of our construction 
that emerges is that the phase space always consists of an A^ x A^ array independently of the 
nature of the group G. We have explored the possibilities for defining Wigner distributions 
on such phase spaces and have shown that this turns out to be possible only when A^ is odd. 
Further, we have shown that for such dimensions, there exist two equally good choices for 
the Wigner distribution and have discussed the relation between them in detail. We have 
also shown that the two Wigner descriptions can be put together into an extended Wigner 
distribution from which they can be recovered by partial tracing. Such a synthesis, however, 
can be achieved only at the expense of introducing overcompleteness into the description. 

The most curious feature that emerges from this work is the importance of being odd. 
We hope to return to this and related questions elsewhere. 
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